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ABSTRACT 

For the very first time, we report the high frequency analysis of Cyg X-1 up to hard X-ray using SPI on-board 
INTEGRAL. After analyzing the possible contribution from the background, and using INTEGRAL archive from 
March 2005 to May 2008, Power Density Spectra (PDS) were obtained until 130 keV. First, we show that their 
overall shape is very similar to that observed at lower energies, as they are well described by sets of Lorentzians. 
The strength of this fast variability (up to 40 Hz) does not drop at high energy since we show that it remains at 
^ 25% rms, even in the highest energy bands. Second, the hard X-ray variability patterns of Cyg X-1 are state 
dependent: the softer the spectrum (or the lower the hardness ratio), the lower the total fractional variability 
and the higher the typical frequencies observed. The strength of total variability as a function of energy and 
state is then investigated. By comparison with simultaneous and published RXTE/PCA data, we showed that in 
the hard state, it remains quite constant in the 2-130 keV energy range. In our softer state, it is also flat until 50 
keV and may increase at higher energy. The implications of this behavior on the models are then discussed. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: binaries — X-rays: individual: 

Cyg X-1 



1. INTRODUCTION 

Black Hole Binaries (BHB), either transient (such as GX 
339-4) or persistent (Cyg X-1), are very variable on long 
time scales (i.e., from weeks to months), in terms of flux, but 
also spectrally speaking: since the first spectra obtained (i.e., 
Tananbaum et al. 1972 for Cyg X-1 and Coe et al. 1976 for 
A0620-00), transitions between two main states were iden- 
tified. The soft state, which spectrum is peaking in the soft 
X-ray band (i.e., ^ 1 keV), is believed to be dominated by the 
emission of the optically thick accretion disk. For that latter 
reason, the soft state is also called "Thermally Dominated" 
(McClintock & Remillard 2006). In the hard state, the Spec- 
tral Energy Distribution (SED) is on the contrary peaking at 
few tens of keV, and is thought to be emitted by an optically 
thin component, the so-called "corona". Its origin, geometry, 
together with the heating processes of this component still re- 
main unclear 

The observed radio behavior of those objects is also de- 
pendent on the spectral state (Homan & Belloni 2005): the 
canonical hard state is usually associated with the presence of 
a compact or sometimes extended steady jet (see e.g. Stirling 
et al. 2001 for Cyg X-1), whereas the observation of relativis- 
tic blobs motion is associated with the intermediate state. In 
the soft state, on the contrary, the radio emission is quenched 
(Fender etal. 1999). 

In order to explain such transitions, a paradigm has been 
set by Esin et al. (1997), where in the soft state the disc is 
extending down to the vicinity of the black hole, whereas it 
is recessed in the hard state at hundreds of gravitational radii. 
The inner part is in turn replaced by an advection dominated 
flow. Observationally, the question to know if the disc is ac- 
tually receding in the hard state is still under debate (see e.g. 
Miller et al. 2006; Reis et al. 2010, but also GierUhski et al. 
2008; Cabanac et al. 2009; Tomsick et al. 2009; Done & Diaz 
Trigo2010). 

BHBs also exhibit strong signatures in the time domain at 
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high frequencies. Their X-ray lightcurves are indeed charac- 
terized by erratic behavior and when analyzed in the Fourier 
domain, their Power Density Spectra (hereafter PDS) can 
show broad band aperiodic noise and/or broad incoherent 
peaks, the so called Quasi-Periodic Oscillations (QPO). The 
strength of those components, together with their typical fre- 
quency evolution, are closely dependent on the state (Homan 
& Belloni 2005). For a complete review of the X-ray binaries 
timing characteristics, see van der Klis (2004). 

Since its discovery in 1964 by an Aerobee rocket, and 
thanks to its persistent brightness due to the wind fed ac- 
cretion process, Cyg X-1 has been extensively studied. The 
aforementioned spectral states are indeed well identified: Mc- 
Connell et al. (2002) e.g., showed that the hard state (which 
is observed at a lower flux in the soft X-ray band), is well 
characterized by the power law component cutting off at 
^100 keV, whereas the soft state (observed at higher flux 
in the low energy bands) is dominated by the disk emission. 
In addition, in both soft and hard states, a high energy tail is 
also observed, which is usually interpreted as the emission of 
non thermal particles (see e.g., Malzac & Belmont 2009 for a 
full SED modeling in presence of magnetic field, and the ex- 
planation of source state transition in terms of particle energy 
distribution variations). 

However, as compared to the other known X-ray binaries, 
Cyg X-1 is in a number of aspect quite different. No relativis- 
tic motion has been observed in the intermediate state (Rush- 
ton et al. 201 1) and the radio emission is not quenched in the 
soft state (Rushton et al. 2011; Zdziarski et al. 201 1). Its soft 
state is peculiar and sometimes called intermediate (Belloni 
et al. 1996). These differences may rely on the high mass 
nature of Cyg X-1 system, where wind fed accretion takes 
place, as compared to the vast majority of other BHBs which 
host low mass companions. 

Cyg X-1 fast variability features were well studied by e.g. 
Pottschmidt et al. (2003), making extensive use of the huge 
RXTE data archive. They confirmed that the PDS could be 
well fitted by sets of Lorentzians, and they showed that the 
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relative contributions of those sub-components were evolv- 
ing with the states. The lower frequency Lorentzian is indeed 
suppressed relative to the second and third Lorentzian during 
the state transitions. Its energy behavior, or rms-spectrum, is 
rather constant in the hard state at ~ 25 — 40% rms. 

There is however an important question which remains 
to address: are these high frequency variability features 
still present in the hard X-rays? At higher energies, a 
strong LFQPO (Low Frequency QPO) was indeed detected by 
SIGMA on board GRANAT (Vikhlinin et al. 1994) in March 
1990 and March 1991. Its typical peak frequency (between 
0.04 and 0.07 Hz) was found to be similar to the broad band 
noise break frequency that was also detected in the power 
spectrum. However its behavior was very variable as the 
QPO, together with the broad band noise had disappeared one 
year later 

The energy dependency of the variability was also exten- 
sively studied in other sources: Rodriguez et al. (2004a,b) 
found that a cut-off near 25 keV was probably present in 
GRS 1915H-105 and XTE J 1550-564 for the LFQPO rms- 
spectrum. Concerning the long time scales variability (from 
days to weeks), Zdziarski et al. (2002) used the RXTE/ ASM 
and CGRO/BATSE ai'chives to study the variability in Cyg X- 
1 from 1.2 to 300 keV. They showed that this very low fre- 
quency variability was also quite high (i.e. 25% rms) and 
state dependent. Examining Cyg X-1 fast variability as a 
function of state, Gieriinski et al. (2010) (hereafter GZDIO) 
showed that the rms spectrum was almost flat for the hard 
state whereas it is increasing with energy in the soft state. In 
intermediate states, it is decreasing with energy. Similar trend 
is observed in other BHBs (see e.g. Gieriinski & Zdziarski 
2005). A tight and strong correlation was also found between 
the Comptonization amplification factor (^h/^s) and the fre- 
quency of the first Lorentzian /i as ^h/^s « /j . 

The models to explain those variability features are still un- 
der debate, principally due to the various patterns observed 
and their evolution with the state. Its origin, i.e., is it coming 
from the disk, the corona, from a transition region or from the 
jet is even not well established. The most probable hypothe- 
sis is that it relies in an interaction between those components. 
Some successful attempts have however been built in order to 
reproduce some of the features observed: Concerning e.g. the 
LFQPO origin, the Accretion-Ejection Instability (Tagger & 
Pellat 1999) was used to explain the LFQPO frequencies and 
the disk geometry evolution (Varniere et al. 2002). Cabanac 
et al. (2010) built a model of an oscillating corona which ad- 
dresses the QPO and the broad band noise typical frequen- 
cies correlation, the modeling of the PDS by radiative transfer 
simulation as well as the qualitative evolution of the frequen- 
cies with states : the softer the SED, the higher the frequen- 
cies. The observed increase of the absolute variability with 
the source flux, so called "RMS-flux correlation" (Uttley et al. 
2005), was finally well modeled by Arevalo & Uttley (2006), 
in the context of fluctuating-accretion model. 

The purpose of this article is to study the fast timing behav- 
ior of Cyg X- 1 at higher energies by using the high sensitivity 
of SPI (Spectrometer Onboard Integral, Vedrenne et al. 2003; 
Roques et al. 2003), as its average effective area is ~ 500 cm^ 
in the 20 keV - 8 MeV band. Its timing ability is indeed un- 
precedented in those high energy bands, and it allowed, e.g., 
absolute timing of the Crab and determination of the hard X-- 
rays to radio delay at ~ 15/is resolution (Molkov et al. 2010). 



2. DATA ANALYSIS PRINCIPLES 

2.1. Event extraction 

A data analysis pipeline was build for timing analysis pur- 
pose of SPI data. This set of Perl scripts is designed to first 
extract the event files, and second the energy spectra asso- 
ciated with given sources in the FOV. The source spectral ex- 
traction is needed for evaluation of the source and background 
contribution to the variability, see hereafter For that purpose, 
we used our local deconv algorithm in our pipeline (see e.g. 
a description of the principles in section 2.3 of Jourdain & 
Roques 2009). The event extraction process is in turn based 
on the algorithm used in Molkov et al. (2010). Some caveats 
described just below were moreover taken into account. 

First, in order to limit contribution from other sources, we 
only retained the pointed observations where Cyg X-1 was in 
the Fully Coded Field Of View (hereafter FCFOV) of SPI, i.e., 
when the instrument axis is closer than 7° from the source. 
Such a choice led us to retain about 1000 SCWs. 

Some SCWs are however affected by telemetry gaps in the 
satellite to earth data transmission. As these lost packets are 
present in less than 10% of the observations, and since it 
affects the Fourier analysis, discarding them was preferred. 
Furthermore, only single events were retained for our study 
(Roques et al. 2003). As stated in Molkov et al. (2010), this 
choice was motivated by the fact that the timing resolution is 
worse with multiple events, together with the difficulty of spa- 
tial information retrieving. Barcentric correction of the event 
files were performed using the standard ISDC barycent 
tool. 

Finally, two energy bands relevant for our study were se- 
lected: the lowest energies were probed between 27 and 49 
keV, whereas the high one is a gather of the 69-90 and 96- 
130 keV bands (abbreviated as 69-130 keV hereafter). Such 
choice is motivated for increasing the source to background 
ratio in the event files by avoiding instrumental lines (see 
e.g. Weidenspointner et al. 2003 for the identification and the 
modeling of those lines). 

Note that the previous caveats are automatically taken into 
account by our SPI timing pipeline, and unless explicitly men- 
tioned, every subsequent errors tabulated or in the text are at 
90%. By default, the plotted ones are at la. The latter confi- 
dence levels for a given value in this article will be expressed 

as Vfl/Me„„„,,„,„^. 

2.2. Evaluation of the noise level in the PDS 

Before calculating any power spectra for a given source, 
it is necessary to check what is the shape and the level of 
the PDSs in absence of the source. For background domi- 
nated instruments, the count rate is high, even when observ- 
ing an empty field region. Hence it is possible to evaluate 
the noise level with good statistics if empty field observa- 
tions are used. We therefore used 45 SCWs from rev. #681, 
corresponding to a total of ~ 156 ks, when INTEGRAL was 
pointing towards mid galactic latitude (RA=19''42'44.23", 
Dec=— 15''04'03.4") , and no significant solar activity was ob- 
served. In order to know how close the noise behaves like 
Poisson statistics, an average empty field PDS (hereafter EF- 
PDS) of rev # 681 was calculated for both 27-49 keV and 
69-130 keV bands. The obtained EFPDS was normalized to 
Leahy (Leahy et al. 1983), such that the comparison with the- 
oretical behavior for Poisson statistics is easier. In that lat- 
ter case, the obtained PDS must be flat when averaged-up on 
infinite number of time frames and a value of 2 is expected 
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Figure 1. Power spectra of the empty field observed during revolution #681 
in the 27-49 keV band (upper panel) and the 69-130 keV band (lower panel). 
Note that except for few points, the PDSs are well fit by a constant (solid 
horizontal line, see text for the parameter of the fits) . 

for purely Poisson noise in Leahy normalization (Leahy et al. 
1983). 

For the 27-49 keV band (see upper panel of Fig. 1) an aver- 
age value consistent with 2 was obtained. More precisely, the 
result of the fit indeed gives Pff_ 27-49 keV = ^■^^^'^'^Imm- 
The statistic for this fit is x^/^ = 1.56 for 44 dof, whereas it 
is x^/'^ = 1-88 when the noise level is set to the theoretical 
value of 2. The moderately good x~ obtained after fitting is 
due to the presence of single peaks in frequency bins, but not 
to an overall trend (e.g. deviation in slope), as either low or 
high frequency bins contribute equally to the x^- Hence, it is 
possible that the peaks observed for the 27-49 keV are only 
due to statistic fluctuations (since our empty field data sample 
is not infinite), but we cannot exclude either instrumental, or 
more probably environmental effects. 

Concerning the 69-130 keV band, it is still consistent with 
a constant close to 2 (see lower panel of Fig. 1), but again the 
fitted value is sUghtly higher as Pp^, 69-i30 keV = 2.00302oq26- 
The statistics is however quite good as x^/i^ = 1 -03 for 44 dof, 
with no significant deviations to the constant obtained. When 
the noise level is set to 2, the statistics give x~/'^ ~ 6.16. 

We tentitatively looked for the reasons of such deviations to 
the pure Poisson noise process, the main one being the pos- 
sible effect of deadtimes on the obtained PDSs. Modeling 



of such effects were performed by e.g. Zhang et al. (1995) 
and used with success on the fine noise level evaluation of 
RXTE/PCA PDSs (see e.g. Morgan et al. 1997). In practice, 
we used the formula given in Nowak et al. (1999) and applied 
it to SPI detector typical values. We obtained a theoretical 
Poisson noise level spanning between 1.98924 and 1.99017 
in Leahy normalization. This theoretical level is therefore 
very different from the value obtained when computing the 
actual EFPDS, and it is possible that some other instrumental 
effects than deadtimes are affecting the EFPDS. Another pos- 
sibility of this excess could be the presence of weak sources 
in the "empty" FOV chosen, which may increase the theoreti- 
cal noise level. Since such deviation from the theoretical noise 
level may also affect the subsequent source PDS calculation 
and its rms evaluation, we built a practical method to over- 
come this effect: even with instrumental effects indeed, the 
source PDS always converge to the noise level at sufficient 
high frequencies. Hence, we determined the empirical noise 
level in each PDS by setting an arbitrary frequency threshold 
/,„,„, and then fit the PDS over / > /,„„, with a constant C„oise, 
which gives finally our practical noise level. Naturally, the 
noise level obtained may contain some source intrinsic vari- 
ability as /„„„ 7^ 00. Hence, some high frequency components 
can be missed or artificially added with that method. We note 
that a similar method (i.e., fit of the noise level at high fre- 
quencies) was used with success on RXTE data for the search 
of twin kHz QPO in neutron stars (Boutelier et al. 2010), and 
we examined the possible effects on an example below (in 
section 3.1). 

2.3. Backgmund dominated PDS 

In the case of background dominated data, corrections have 
to be applied in order to evaluate the intrinsic power spec- 
trum of the source. We indeed used the raw count rate de- 
tector lightcurve as input for calculating the raw power spec- 
trum Plot, Leahy It is then easy to demonstrate that the intrinsic 
power spectrum Psrc,rws expressed in "rms" or "Myiamoto" 
normalization (Belloni & Hasinger 1990, or Miyamoto et al. 
1991) follows: 



[P' 



tot.Leahv ' 



' ^noisej 



(B + s) 



(1) 



where (5') is the average source count rate, and {B + S) the av- 
erage detector count rate (see also Vaughan & Nowak 1997). 
For coded mask instruments, the value of (5) is obtained by 
using the usual deconvolution algorithms for spectral extrac- 
tion, and the average detector count rate {B + S) is directly 
given by the event files. 

3. RESULTS 

For the power spectra computation, we arbitrarily set each 
segment duration to AT = 647 s, fixed the binning time such 
that the highest frequency available is /Aiy^y = 810 Hz and re- 
binned logarithmmically the power spectra obtained. For that 
purpose, we used the standard powspec XRONOS tool to 
compute the mean power spectrum {Ptot. Leahy)- We then used 
ISIS' version 1.6.1-26 (Houck & Denicola 2000), as well 
as routines publicly available from the S-lang/ISIS Timing 
Analysis (SI TAR)- for PDSs post-processing (fitting, time 
lags and rms evaluation). The noise level was finally eval- 

' See http://space.mit.edu/CXC/isis/index.html 
-http://space.mit.edu/CXC/analysis/SITAR/ 
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Figure 2. Power spectrum of Cyg X-1 calculated from rev #301 to 684 (ob- 
servations where EXO 2030+375 was flaring were removed) in the 27-49 
keV band. The intrinsic average fractional variability is rms ~ 29 ± 3%. 4 
Lorentzians were used to fit the spectrum. See first column of table 1 for the 
parameters of the fit. 

uated as prescribed in the previous section, i.e., we fitted 
the high frequency tail of each PDS by a constant. The 
minimum frequency chosen for the fit of Cnoise was set to 
fmin. noise = 100 Hz. We also rcbinncd again the obtained 
PDSs at high frequencies for better statistics which explains 
the non equal bin spacings in the PDS plotted. 

3.1. State averaged power spectra ofCygnus X-1 

We first focused in our study on state-averaged power spec- 
tra computation. The following PDS hence mixed different 
physical state of the source, but this choice was motivated 
primarily for increasing the statistic at high frequencies in the 
highest energy band probed. We also found that averaging up 
the observation from revolutions #301 to 684 were the best 
trade off for increasing the statistic in the 69-130 keV bands. 
Earlier observations (especially revolutions #79 and 80) are 
indeed quite soft in the SPI energy range, and therefore de- 
creases the signal to background ratio of the source in that 
band. 

Noise subtracted power spectra of Cyg X-1 in "rms" nor- 
malization are plotted in Fig. 2 for the 27-49 keV band and 
in Fig. 3 for the 69-130 keV band. They were calculated 
from rev #301 to 684, but during this period the accreting pul- 
sar EXO 2030H-375 located 7.2° from Cyg X-1 was flaring. 
Hence, as the total SPI FOV is 35° (16° for the fully coded 
FOV), and to prevent any contribution in the calculated PDS 
from that source, we removed from our data sample the ob- 
servations corresponding to the dates of the flare, i.e from rev 
#452 to 486. Disentangling any contribution from this source 
in the PDS is the purpose of the appendix. 

Our first result is that the overall shape of the noise sub- 
tracted PDSs in the 27-49 keV band shown on Fig. 2 looks 
very similar to the ones obtained with softer X-ray monitor- 
ing of the source by e.g. RXTE (see e.g. Lin et al. 2000; 
Nowak et al. 1999). Fitting the obtained PDS by a constant 
gives in turn a very bad reduced x~ since x~/^ = 37.8 for 44 
dof. We therefore chose to fit the PDSs as it is usually done in 
the literature, i.e., with sets of Lorentzians L, (Nowak 2000; 
Befloni et al. 2002; Pottschmidt et al. 2003), i.e.; 



LiH 



1 



2Rf 




Frequency (Hz) 

Figure 3. Power spectrum of Cyg X-1 calculated from rev #301 to 684 (ob- 
servations where EXO 2030+375 was flaring were removed) in the 69-130 
keV band. (Up): 3 zero centered Lorentzians were used to fit the PDS. 
(Down): 2 zero centered and 2 "QPOs". See second column of table 1 for the 
parameters of the fit. 

with /, the resonance frequency of the Lorentzian, Aft the 
full frequency width at half maximum value, and R/ a normal- 
ization constant. We then define the quality factor Qi of this 
Lorentzian as; 

the normalization Rj is linked to the rms of this Lorentzian by 
integrating L, over the whole frequency domain, which gives 
(see also Pottschmidt et al. 2003): 



'1 atan{2Qi) 
rms, — Ri\l — - 



(4) 



ttA/,- 



1+4 



m 



(2) 



In the following text, fitting the PDS with "QPO" refers to 
Lorentzian with /, 7^ 0. However, for sake of simplicity, we 
generally used by default zero frequency Lorentzians, i.e., 
fi = 0. This choice was motivated for decreasing the num- 
ber of free parameters in the fit while giving reasonable x^ 
values. In that case, Qi = and the only relevant parameters 
are A/, and the normalization R, ~ rms; \/2. 

For example, fitting the 27-49 keV PDS with only three 
Lorentzians gives a reduced x^ /v = 1.53 for 34 degrees of 
freedom (dof), leaving a high frequency tail unfitted. The 
addition of a fourth zero centered Lorentzian gives in turn 
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Table 1 

Revs #301 to 684: results of the PDS fit by 3 or 4 zero 
centered Lorentzians. 



Energy band 


27-49 keV 


69-130 keV 






(a) 


(b) 


A/o (Hz) 


n nn94°°°''3 

U.UU/4g Q(,Qg7 


0.0038(, QQQ, 


0.0024 (/) 


rmso 


U-U3^0.026 


0.031°;°^3 


r. 0770035 
U-^'-^ '0.016 


A/i(Hz) 




0.25»:J5 


0.15 (/) 


rms ] 


n i3''i'» 
?-o§:ii 

^•^^0.85 


0.13°:ii 


o-iii!:lo 


A/2 (Hz) 




1-39 (/) 


rms2 


U.ligi, 




o.o9°:» 


A/3 (Hz) 


38 1 
23^^^ 
0.29!];i 


32^« 


38.1 (/) 


rms^ 


0.26°:1° 

o.3oi|:i 


n 76° 32 


rmstoi/mean 


xV^ 


0.87 


1.04 


1.08 



Note. — (/) stands for a parameter value fixed in the fit, (p) for a 
parameter that pegged at a limit during the error calculation. 

'^' Values obtained when all parameters of the fit are let free to vary. 
' ' Values obtained when tying the Lorentzian frequencies to the value 
obtained in the 27-49 keV band. 

X^/i/ = 0.72 for 32 dof. It corresponds to a F-test Null Hy- 
pothesis Probability (hereafter NHP) of 2.4 10"^ for this ad- 
dition (see however Protassov et al. 2002 for the limited uses 
of the F-test), which gives good confidence in the presence 
of such fourth component. However, since the presence or 
the absence of this fourth component may be affected by an 
under estimation of the noise level, we also reanalyzed the 
data in a worse case: we set the minimum frequency for the 
noise level fit to /„,„, „oise — 0.005 Hz instead of 100 Hz. With 
that choice, the noise level obtained will always increase, 
as it may include some source contribution as well. There- 
fore, after subtracting this new noise level, we re-fitted the 
PDS with three or four Lorentzians. The obtained x^ values 
are naturally worse, as compared to the previous case where 
fmin, noise = 100 Hz, as they are dominated by systematic neg- 
ative residuals at high frequencies: we obtain x^/'^ = 1-83 
with 34 dof for 3 Lorentzians and x'/i^ — 1-29 with 32 dof 
for 4 Lorentzians. Hence, even with this extreme evaluation 
of the noise level, the F-test NHP for the addition of the fourth 
Lorentzians is equal to 1 .5 10^^. This gives in turn good con- 
fidence for the presence of such high frequency tail in our data 
set. 

The fit parameter values, together with their corresponding 
errors are displayed in Table 1 . The errors on each fitted pa- 
rameter were calculated leaving all the others free to vary, un- 
less mentioned explicitly. The four fitted Lorentzians frequen- 
cies values span from 2.10^^ to 35 Hz and are logarithmically 
spaced by a factor ~ 5 to 20. Besides the very low frequency 
component of frequency A/o which account for less than 3% 
of the variability, the low frequency Lorentzians (Li and L2) 
contribute to the total rms by the same order of magnitude 
amplitude, i.e., 12-13%, and it increases to 23% for the high 
frequency (L3) Lorentzian. The overall Cyg X-1 variability in 
the 27-49 keV band is therefore estimated to 29.4 ± 2.8% rms 
from revs #301 to 684. 

For the 69-130 keV band, the PDS qualitative shape was 
showing even more deviation to the powerlaw shape with 
clear defined "humps" (see Fig. 3). Again, the PDS fitted with 
a simple constant gives x^/i' ~ 5.26 for 44 dof. When fitted 
with a model of 3 zero-centered Lorentzians, a quite good re- 
duced x^/i' ~ 1-04 for 31 dof is obtained. The addition of a 
fourth Lorentzian is not required as it gives roughly the same 



X^, and the corresponding F-test NHP is in turn 38%. Note 
that the Lorentzian which is not required as compared to the 
27-49 keV is the "mid" frequency ones (with frequency A/2). 
Since the errors on the parameters are bigger than the ones 
obtained for the 27-49 keV, we also tried to fit the 69-130 keV 
PDS by tying the typical frequencies of the Lorentzians to the 
values obtained at low energies. This in turn gave a similar x^ 
value ix^ jv — 1.08 for 34 dof), and the total rms values are 
equivalent, typically around 30% rms. It is also interesting 
to note that when the frequencies are tied in the 69-130 keV 
fit, the rms of each sub-component is consistent with the one 
found for the 27-49 keV band. 

As the the 69-130 keV PDS exhibits humps and espe- 
cially at low frequencies, fitting them with non zero centered 
Lorentzians (i.e., QPOs) was also tried (see lower panel of 
Fig. 3). The reduced x^ value can be reduced down to 0.80 for 
17 dof by this way, if a model composed by 2 zero centered 
Lorentzians and 2 "QPOs" is used. This alternative model 
does not affect the overall variability estimation as in that 
case rmsioilmean = 0.28q24. The values of the quality fac- 
tor Qi are Qx = 1.24™° and Q2 = 0.62^;4j for both "QPOs" 
fitted. Hence, since Q <2, those humps found in the 69-130 
keV PDS may be qualified as "peaked noise" (see however 
e.g. Belloni et al. 2002 and reference therein for a criticism 
about such a strict QPO/peaked-noise boundary). The fre- 
quencies of both peaked noise are in turn f\ = 0.033o 02? Hz 
and/2=0.22i]:25Hz. 

The results of the 69-130 keV PDS zero centered 
Lorentzians fits are given in Table 1. In summary, the PDSs 
shows strong variability from lO^"* to at least a few tens of Hz, 
reaching a level of roughly 30% rms, even at those high ener- 
gies. The low frequency component (Li or L2 when present) 
contributes for 10 — 15% to the rms. In contrast, the high fre- 
quency tail contributes for at least 26%. 

3.2. Evolution of the rms 

We estimated our overall rms by integrating each modelled 
PDS. In this section, we studied the evolution of the variability 
and its dependence with the state or the energy. 

3.2.1. RMS as a function of time 

The evolution of the variability as a function of time was 
first examined. In order to obtain sufficient statistics for eval- 
uating the rms, but also its evolution with time, a one rev- 
olution typical time scale was found to be the best trade-off 
for grouping the data together. For the same reasons, we also 
excluded the revolutions where Cyg X-1 was in the FCFOV 
during less than 20 SCWs. Hence, 16 different revolutions 
were selected, spanning from revs #79 (June 2003) to 684 
(May 2008), and corresponding to a total of 2.2 Ms of data. 
Noise subtracted PDSs were then computed for each revolu- 
tion in the same way as done in section 3.1. Only 1 or 2 zero 
centered Lorentzians were necessary to fit each of them well, 
given the lower statistics compared to the previous grouped 
analysis. In order to compare the timing parameters evolu- 
tion with the spectral ones, we also built deconvolved spectra 
for each selected revolution. Using those spectra, the source 
count rates in the 27-49 keV and the 69-130 keV bands were 
finally computed. In order to get a quick look on the spectral 
shape of the source, we also defined a SPI hardness ratio as: 
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Table 2 

Results of the correlation tests between the fractional rms in the 

27-49 keV band and count rates or HRs. Prob. refers to the null 

hypothesis probability under the considered test (Spearman's p 

or Kendall's r rank coefficients). 
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SPI count rate 
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Figure 4. Evolution of Cyg X-1 timing and spectral parameters from our sample. Left panel: Evolution of rmstot/mean in 27-49 keV band with time. Center 
panel: SPI hardness ratio (HR) as a function of the derived ASM photon index Tasm ■ Right panel: rms vs SPI HR. Straight fines sfiow the fits of the points in 
the middle and right plots by a linear function. 90% errors plotted. 

3.3. Variability in Hard State vs Soft State 

In the following section we examined the fractional vari- 
ability of Cyg X-1 as a function of state. For defining those 
states, we used two different methods. The first is based on 
the photon index Tasm obtained in soft X-rays, and the sec- 
ond one is based on the hardness ratio obtained in the hard 
X-ray band {HRspi). Both are anti-correlated as demonstrated 
in the above paragraph. However, due to the scatter of Tasm 
vs HRspi plot, grouping of the observation based on one or 
the other criterion lead to different samples, and hence, as we 
will demonstrate, different results. 

3.3.1. States defined by the ASM 

The first state definition is based on the derived value of 

Tasm by: 

• Tasm < 2 corresponds to a hard state (HSasm), whereas, 

• Tasm > 2 corresponds to a soft state (SSasm)- 

The same process as described in section 3.1 was applied to 
the grouped HSasm and SSasm event files. The 27-49 keV 
PDSs are plotted in Fig. 5. 

For the HSasm, the 27-49 keV PDS is remarkably well 
fit by three Lorentzians ix^/v ~ 0.29, whereas x^/'^ = 1-66 
with only 2 Lorentzians). The typical frequency values are 
A/i = 0.26[!i^ H2 jjnd ^y^ ^ 5 39.5 ^^ g^^j ^^ ^Q^^l fj.^^._ 

tional variability is 18^^%. The fit values obtained in the 69- 
130 keV band are similar though the typical frequency found 
for A/i = 0.44q 25 Hz is a bit higher and the overall variability 
decrease to 12jy%. 

For the SSasm, only two Lorentzians are needed to fit well 
(x^/v = 1 .66) the 27-49 keV PDS. The overall variabiHty is at 
11 ± 2%, hence lower than in HSasm, but the first Lorentzian 
frequency increased as A/i = 0.50q 27 Hz. In the 69-130 keV 
band, the variability remains at a level of 1 1 ± 5%. 



The rms in the 27-49 keV band as a function of time is shown 
on Fig. 4. The overall trend is that the total fractional rms 
evolves with time from about 10% to 35%. This variable be- 
havior led us to examine any possible correlation between the 
calculated rms and the spectral parameters. 

3.2.2. RMS - hardness ratio correlations 

The fractional rms as a function of SPI hardness ratio is 
plotted in Fig. 4. The overall clear trend is that the frac- 
tional variability tends to correlate with HRspi: rank coeffi- 
cients were computed (see Table 2) and gave a NHP equal 
2.10^^ (2.10^"* for the Kendall's test), hence the rms is well 
correlated with HRspi- x" fitting of the previous relations 
were also tried, in order to take errors into account. Sim- 
ple linear functions were used for those fits (straight lines in 
Fig 4) and the correlation between the rms and HRspi is still 
confirmed by this way. The line slope is indeed S,ms-HRspi = 
0.64 ±0.11 rms. 

3.2.3. ASM photon index vs SPI liardness ratio 

The important question of how HRspi relates with the state 
definition usually based on soft X-ray criteria was then ex- 
amined. For that purpose, we retrieved from the RXTE- ASM 
archive^ the simultaneous count rates of the source in the low 
(A, 1.5-3 keV), medium (B, 3-5 keV) and high (C, 5-12 keV) 
standard energy bands. Those count rates were then converted 
into fluxes and 3-12 keV photon indices (Tasm) of the source, 
following the method described in Zdziarski et al. (2002). Fi- 
nally, HRspi was plotted against the obtained Tasm (see cen- 
ter Fig. 4). It appears that HRspi is anti-correlated with Tasm 
(Spearman's rank coefficient is here p = —0.49), as when fit- 
ted with a straight line, the HRspi vs Tasm gives a negative 
slope of Shrspi vs Tasm = ^ 1 -95 12 72- The y axis intercept value 
for this straight line is in turn Thrsp,=() = 3.272 ^y 



3.3.2. States defined by SPI 

Since the errors on Tasm are not negligible (see center panel 
of Fig. 4) and make the previous observation grouping ques- 
tionable, we also tried using the SPI hardness ratio for defin- 
ing two different states: 

• HRspi > 0.75 corresponds to a hard state (HSspi), 
whereas. 



' http://xte.mit.edu/asmlc/srcs/cygx 1 .html 



HRspi < 0.75 corresponds to a soft state (SSspi)- 
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Figure 5. Cyg X-1 PDSs calculated when the source was hard (left panel) and soft (right panel) according to the ASM criteiion in the 27-49 keV band. 
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Figure 6. Cyg X-1 PDSs calculated when the source was hard (left panel) and soft (right panel) according to SPI hardness ratio criterion in the 27-49 keV band. 



This value of 0.75 was chosen because of the clear separation 
visible in HRspi between two sets of data. Such a separation 
would correspond to a boundary value in ASM photon index 
of 1.80 according to the previous HRspi vs Tasm fit, but is 
compatible with the previous boundary of Tasm = 2, when 
considering the typical errors on this fit (ArH«s,>,=o = 0.52 at 
90%). 

The 27-49 keV band PDSs for both SPI defined states are 
plotted in Fig. 6. For the 27-49 keV PDSs, 3 Lorentzians 
were necessary to fit the HSspi PDS, whereas only two where 
used for the SSspi one. The addition of a second component 
was tested for the HSspi and the corresponding F-test NHP 
is equal to 6.3 10^^. The addition of a third one gives a F- 
test NHP equal to 1.5 10^1 For the SSspi, the addition of a 
second component (with parameter A/i) is highly required as 
it gives a F-test NHP equal to 2.7 10^^. The addition of a 
third Lorentzian in the SSspi PDS would is unlikely since the 
corresponding calculated F-test NHP is equal to 8%. For the 
69-130 keV, only two Lorentzians were used for HSspi and 
one for SSspi- 

The results of the fits of both states, are given in Tables 3 
and 4. It shows first that the HSspi overall variability is quite 
high in the 27-49 keV band (26 ± 3%), and dominated by a 
high frequency (~ 20 Hz) component. As this latter com- 
ponent is not needed in the 69-130 keV PDS, the variability 



Table 3 

SPI Hard state: results of the PDS fit by 3 Lorentzians. 
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Note. — (/) stands for a parameter value fixed in the fit. 
^^' Values obtained when all parameters of the fit are let free to vary. 
^ ' Values obtained when tying the Lorentzian frequencies to the 
value obtained in the 27-49 keV band. 

drops down to 13 ±2%. Tying the frequency values in the 
69-130 keV PDS to the one obtained in the 27-49 keV band 
lead however to a rms value of 24;^% consistent with the one 
obtained at low energy. 

In the SSspi, the frequency of the first Lorentzian increases 
up to A/i = 1.3 Hz. The variability is lower than in HSspi 
(~ 15%) in the 27-49 keV band, but it increases in the 69-130 
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Table 4 

SPI Soft state: results of the PDS fit by 1 or 2 
Lorentzians 



Table 5 

Comparison between SPI and PCA PDSs: results 
of the fit 
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keV band (up to ~ 20% or even ~ 30% when the frequencies 
are tied to the ones obtained at low energy). 

Comparing to the previous state grouping based on ASM 
criterion, it is interesting to note that the typical frequencies 
obtained here are more different between the hard and the soft 
state. For example the frequency A/i in SSasm is only twice 
the one in HSasm, whereas there is a factor 7 for A/i between 
SSspi and HSspi- A difference in the values of A/2 is also 
noted between HSspi (where A/2 ~ 20 Hz is high) and HSasm 
(where A/2 ~ 5 Hz is low). 

3.4. Comparison with RXTE/PCA power spectra 

We also compared our INTEGRAUSPl PDSs with 
RXTE/PCA ones. The only contemporary RXTE observa- 
tions found for our dataset spans from revolution # 628 to 
630 for INTEGRAL. During those observations, the average 
ASM photon index is hard as Tasm = 1.8 ±0.15. We re- 
trieved and analyzed in turn the whole observation #93121 
for RXTE/PCA. The PCA observation is simultaneous to IN- 
TEGRAL ones, but the duration of the latter is larger for in- 
creasing the statistics. For the PCA observations, event files 
were processed in the standard way with seextrct and 
background lightcurves were calculated with pcabackest 
and saextrct. We finally computed the PCA PDS us- 
ing powspec on the 20-40 keV background substracted 
lightcurves. 

Both PDSs are plotted on Fig. 7. It is first clear that 
the overall shape is similar from one instrument to the other. 
When fitted with two Lorentzians, the typical frequencies are 
indeed compatible, together with their normalization (see Ta- 
ble 5). The only noticeable difference is the presence of a 
very low frequency complex present in SPI PDS whereas it 
is absent in the PCA ones. This absence could reflect the 
shorter duration of the PCA observation, and there would be 
not enough variability at frequencies comparable to its dura- 
tion. In any case, setting the SPI PDS fit parameters to the 
value obtained for the PCA PDS still gives an acceptable fit 
as in that case x^/i' = 1.35. 

Finally, for the 69-130 keV band, due the lower statistics, 
only one Lorentzian is required for fitting the PDS {y^ jv = 
1.04). The typical frequency is consistent with the ones ob- 
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and 



tained at lower energy since A/j 59- 
its variability is rms\, 69-130 keV = 0.12qo9. As the lower en- 
ergy PCA PDSs were also analysed for this observation, for 
the first time a 2-130 keV high frequency rms spectrum of 
the source was build. It is plotted in Fig. 8. Note that 
HRspi = 0.84±0.01 for this observation, hence Cyg X-1 was 
in a hard state at this time. The rms spectrum is flat until 
49 keV (- 20%), and decreases in the 69-130 keV bin. Fit- 
ting the rms spectrum with a nearly flat powerlaw of index 
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a = 0.015qoq2 gives a bad yi^ jv = 2.4, and the addition of a 
high energy cutoff improves the fit as 'y^ jv ~ 0.95. The ftest 
NHP for such addition is however 1.5%, hence its presence 
is rather marginal, together with the lack of constraints on the 
fitted folded Efoid = 31^p keV and cutoff £„„ = 84{f keV 
energy values. The latter statistic is also dominated by the 
last bin (69-130 keV rms), as the 2-49 keV rms spectrum is 
well fit by a simple powerlaw {jc /v = 0.94). 

3.5. Time lags 

We looked for time lags between the low energy band (27- 
49 keV) and the hard one (69-130 keV). For that purpose, 
we used the sitar_lags function provided by the SI TAR 
package. Lags as a function of frequency were then obtained 
and are plotted in Figs. 9. It appears first that the obtained lags 
are dominated by the Poisson noise due to the background 
contribution. 

In the hard state, a strong rebinning of the cross-power 
spectrum tends to show that the calculated lags are consis- 
tent with or slightly positive (i.e., 69-130 keV band lags the 
27-49 keV band) at low frequency (i.e < 10~^ Hz). The effec- 
tive lag calculated is e.g., Ar = +18±12s(1(T errors) in the 
[1.5,9.2]10~^ Hz frequency band. Between 1 and 10 Hz, an 
upper limit of At < +30 ms can be set. This limit is compat- 
ible with the positive lags observed in softer X-ray in the hard 
state, where it is typically Ar ^ 1 — 4 ms (see e.g. Pottschmidt 
et al. 2003). When integrated over a wide frequency range, 
cross-power spectra are consistent with lags at (90%) for 
this hard state: we indeed obtain Ar = +6.8 ±6.4 s {la er- 
rors) in the [00.0015,1.5] Hz frequency range. In our soft 
state, the calculated lags seem also positive and slightly better 
constrained at low frequencies. For example, Ar = +19±7 s 
in the [0.92,5.5] 10^~ Hz frequency band. As for the hard 
state, an even stronger rebinning of the cross-power spectrum 
results in obtaining a value consistent with lags (at 90%) as 
in that case Ar = 35 ± 29s (Ict errors) in the [00.0015, 1 .5] Hz 
frequency range. Those results are consistent with those ob- 
tained by Torii et al. (2011), as their cross-correlation func- 
tions between the 10-60 keV and 60-200 keV peak at lag. 

4. DISCUSSION & CONCLUSIONS 

SPI on-board INTEGRAL was used successfully to study 
Cygnus X-1 variability at high energies. The high stability of 
the instrument allows to probe Cyg X-1 variability up to 100 
Hz and down to 1 mHz. The first and main result of this study 
concerns the evolution of the variability at high energies: it is 
indeed demonstrated that Cyg X- 1 is varying at a high level 
(-25%) up to 130 keV. 

The shapes of the high energy PDS obtained are similar 
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Figure 7. Cyg X-1 PDS for SPI in the 27-49 keV band (left panel) as compared to contemporary PDS in the 20-40 keV with the PCA (light panel). See text for 
the result of the fit. 
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Figure 8. Cyg X-1 rms spectrum for the contemporary INTEGRALIRXTE 
observations. Crosses = PCA, diamond = SPI data. 90% errors plotted. 

to that obtained at lower energies with other instruments as 
it is well described by subset of up to 4 Lorentzians. High 
frequency variability {^ 40Hz) is observed in the 69-130 state 
averaged PDSs, showing for the first time that fast variability 
is present at high energy. 

The parameters of the Lorentzians used are also evolving 
with the source state. For example, in order to explain the 
differences between the HSspi and SSspi frequencies found, 
the simplest solution is to associate those frequencies with 
given physical components, at certain typical radii in the disk 
or in the corona. In our case, the first Lorentzian is shifted 
from 1 .6Hz in the SSspi to 0.22Hz in the HSspi- If this com- 
ponent is linked to Keplerian motion, it would imply that 
the typical radius for this component increases by a factor 
(0.22/1.6)"^/^ = 3.8 between the soft and hard state. In the 
framework where the variability is driven by waves propagat- 
ing in the corona (see. e.g. Cabanac et al. 2010), the typical 
break frequency scales as the inverse of the external corona 
radius. Hence in that latter framework the corona external ra- 
dius would be 7 times larger in our HSspi than in the SSspi- 
It is interesting to note that the frequency shifts for /i follow 
the difference between the hard and soft state fundamental 
planes found by Kording et al. (2007): this latter is based on 
the fundamental plane found by McHardy et al. (2006) which 
links, for SS BHB and AGN, the mass accretion rate m, the 
first Lorentzian frequency /i and the black hole mass M with 
/i X M ex lit. Moreover, Kording et al. (2007) have shown that 
this relation still holds in the HS, at price of changing the pro- 



portionality factor. Based on an extensive analysis, they cal- 
culated that this factor was equal to 0.3 for HS BHB, whereas 
it was 2.2 in the SS. In our case, we note that the ratio be- 
tween A/i. soft and A/i. hard is equal to 7, which is therefore 
very close to the proportionality factor ratio in both funda- 
mental planes. Such different tracks for HS and SS states in 
the fundamental plane is interpreted by Kording et al. as a 
sign that not only the accretion rate is driving the accretion 
process evolution. An additional parameter must be invoked 
to render this two different tracks in the Hard and Soft states. 

The rms evolution with the hardness of the source that we 
observed at high energy is similar to the rms-hardness corre- 
lation seen in many XRB (see e.g., the Hardness-RMS Dia- 
grams in Belloni 2010). It reflects the known relation that the 
source gets more variable when it gets harder (see however 
some high variable soft state in GZDIO, Fig. 2b). Our results 
support the fact that this relation still holds at high energies, 
giving strict constrains on emission models that intend to re- 
produce the observed variability. 

The rms level we obtained (~ 25%) together with the fre- 
quency of the first Lorentzian (A/i ~ 0.2 Hz) in our hard 
state is compatible with the majority of the hard states found 
at lower energies with the PCA in GZDIO (see e.g. Figs 2b 
and 7a of that paper). If in GZDIO the hardest states rms 
reaches 35-40%, our hard state level, PDS and typical fre- 
quencies are compatible with bit softer ones. In particular, our 
hard state PDS (Fig. 6) is close in shape to the PCA PDS for 
observation 60090-01-19-00 (Fig 6b of GZDIO), for which a 
rms spectrum is computed: the variability is rather constant 
around 25-30% rms on the 2-30 keV band, and adding our 
two high energy point would show that the rms remains at 
roughly the same level in the 2-130 keV band for this hard 
state. This rather constant rms spectrum in the hard state 
is emphasized by our simultaneous PCA/SPI analyse of the 
source, where the source was hard. If a break in the hard 
X-ray (~ 80 keV) is possible in this rms spectrum, its signif- 
icance is rather marginal. It is also possible that the low du- 
ration of this observation and hence the lack of statistics does 
not allow to constrain a putative high frequency tail, which 
could account for extra variability. Independent results from 
Torii et al. (201 1) obtained with Suzaku showed that the hard 
state of Cyg X-1 rms is indeed rather constant in the 12-200 
keV energy range. 

Our SSspi can be similarly associated with PCA observation 
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Figure 9. Lags between the 69-130 keV and the 27-49 keV band as a function of energy for the hard (left panel) and soft state (right panel). Filled diamonds 
stand for 69-130 keV band lagging the 27-49 keV bands, i.e., positive lags, whereas crosses stands for negative ones. 



id 70104-01-03-00 analysed by GZDIO. The rms spectrum 
is in turn also quite constant between 2-35 keV, though its 
level is lower than in the previous hard state. Our SPl level 
of variability in the 27-49 keV band is consistent with the 
one obtained in the highest energy band of this PCA obser- 
vation (15-20%) and it increases in the 69-130keV (20-30%). 
An interesting interpretation of the rms spectra has been at- 
tempted by Gierliriski & Zdziarski (2005) and more recently 
in GZDIO, with the following principles: any lightcurve vari- 
ability measured by its fractional rms in a given energy band 
may be due to the variability of at least one of the parameters 
Pi of the average spectrum. Gierlihski & Zdziarski (2005) and 
GZDIO used the combination of a multicolor disk (diskbb 
Mitsuda et al. 1984), the eqpair hybrid comptonization 
model (Coppi 1999) and a reflection component (Magdziarz 
& Zdziarski 1995), in order to model the average spectrum F. 
The effect on the output rms spectrum of each parameter vari- 
ation is then investigated. Qualitatively, if a parameter /?,■ only 
impacts on the normalization of a model, the spectrum will 
be shifted by a certain factor, independently from the energy. 
Hence, the variability of this parameter will lead to a flat rms 
spectrum. On the contrary, if another parameter pj (e.g. the 
photon index F) changes the slope of the spectrum around a 
pivot, the factor by which the spectrum is shifted will not be 
energy independent: the further the working energy band is 
from the pivot, the higher the rms. 

For the hard state, a flat rms spectrum is observed for XTE 
J1550-564 and XTE J1650-500 (Gierlihski & Zdziarski 2005) 
and Cyg X-1 (GZDIO). In those papers, it is proposed that the 
flat rms spectrum in the hard state may be reproduced by the 
number of soft photon variability or fluctuation of the accre- 
tion rate. Given that our rms spectrum behaviour remains flat 
until 130 keV in the hard state, we demonstrate that this inter- 
pretation still holds at high energies. 

For our soft state examined here, the rms spectrum which 
is flat until 50 keV, may then increase at high energies. Fol- 
lowing (Gierlihski & Zdziarski 2005), variability of the in- 
jection Lorentz factor r,„y may reproduce such rms spectrum 
increase at high energies, but a maximum should also be ob- 
served around 5-10 keV. 

The above explanation of rms spectra however first depends 
on the underlying average model chosen, and second on the 
fact that only one parameter variation is considered. It is 
hence possible that some of the parameters are linked by un- 
derlying physical processes, such as e.g. the luminosity of the 
hard and soft components when m is varying. As an exam- 
ple, Cabanac et al. (2010), in the framework of an oscillating 



corona model, examined the impact of an adiabatic oscillation 
on the rms spectra. Both the corona temperature and the op- 
tical depth are then varying, which in turns leads to a Comp- 
ton y parameter variation. They predicted the presence of a 
minimum in the rms spectrum, due to the pivoting around the 
average spectrum. The position of the pivot depends on the 
relative peak between the soft photon source and the comp- 
tonized component. 

We finally note that the high energy band power spectra 
that we obtained exhibits peaked-noise components, one of 
wich (the lower frequency one) is consistent with the QPO 
detected with SIGMA by Vikhlinin et al. (1994) in terms of 
peak frequency value. Those components are sharper in the 
higher than the lower energy PDSs (comparing e.g. Figs. 2 
and 3). This effect is also observed at lower energies with 
RXTE (see e.g. Fig. 2 of Cui et al. 1997). If confirmed, 
it may therefore indicate that the Lorentzians quality factor 
increases with energy. In other words, the coherence of the 
underlying variability processes would increase with energy, 
giving constraint on the emission models. 
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APPENDIX 
CONTRIBUTION OF OTHER SOURCES IN THE FOV : DECONVOLUTION ALGORITHM FOR PDS WITH SPIROS 

Our deconvolution algorithm is based on the fact that the Leahy Power is a hnear function of the count rate. As a consequence, 
any deconvolution algorithm can be used for the discrimination between the sources PDSs, if the PDSs are Leahy renormalized. 
Methods using either PIF (Partial Illumination Fraction of a detector) or usual deconvolution algorithm can be used. PIF methods 
(see e.g. an example applied to IBIS/ISGRI data, Cabanac et al. 2006) are quite relevant when the detector number of pixels is 
large. For SPI, since the pixel number is low (19), usual deconvolution algorithm was preferred here. Hence, first for each SCW, 
19 different PDSs are calculated (after retrieving and selecting the events according to the energy band chosen). Second SPIROS 
(ISDC standard deconvolution algorithm Skinner & Connell 2003) is used in spectral mode, but the input files are modified: 
counts per energy band and detector are usually employed, but they are here replaced by the Leahy Power per frequency bin 
and detector. The same process was applied to the background files. We can therefore summarize this process by the following 
substitution in the input files; 



Energy bin - 

Detector Counts - 

Back. Counts - 



► Frequency bin 

- Leahy Power x Livetime 

'Noise Leahy Power X Livetime. 



(1) 
(2) 
(3) 



The multiplication of the Leahy Power by the actual live-time of the observation is due to the fact that the Leahy Power scales as 
the count rate, whereas the usual deconvolution algorithm uses the total counts. We finally ran SPIROS in background handling 
mode 3, and obtained the deconvolved PDSs for for Cyg X-1 and EXO 2030H-375. They are plotted in Fig. 10. In our sample, 
the EXO 2030H-375 pulse and its 2 first harmonics are clearly detected in the associated PDS whereas it is absent in Cyg X-1 
PDS as expected. A fit of the pulse and its harmonics by narrow Lorentzians gives a peak frequency of fpuise = 24 mHz, which 
corresponds to a pulse period of TpuUe = 41.6s. This value is consistent with the one derived by other instruments/observatories 
for the 2006 outburst (see e.g. Klochkov et al. 2007). None of these features are present in the Cyg X-1 PDS, which gives good 
confidence in our deconvolution algorithm. 

We also fitted the continuum with one zero centered Lorentzian. During these observations, we estimated the fractional 
variability of Cyg X-1 to be rmsi = rms,ot/mean = 252*1%. ^'^^ EXO 2030H-375, fitting the PDS continuum with a single 
zero centered Lorentzian gives a rather bad fit {x^/v — 1.91), but gives an evaluation of the low frequency variability with 
rms\ ~ 3531%. Finally, the Leahy level obtained for EXO 2030H-375 during its giant outburst is comparable to the one obtained 
for Cyg X-1. Hence, given that EXO 2030H-375 flux is usually ~ 10 times fainter"*, we confirmed that our first approximation of 



e-g- 



the SwiftlQKi lightcurve 



http://heasarc.gsfc.nasa.gov/docs/swift/results/transients/EXO2030p375/. 
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Figure 10. Deconvolved (and noise subtracted) rms normalized PDSs of Cyg X-1 (up) and EXO 2030+375 (bottom) calculated from rev #452 to 486 in the 
27-49 keV band. Note the pulse fundamental and its harmonics well associated with EXO 2030+375, and not present in Cyg X-1 PDS as expected. 

neglecting this source in the raw PDS calculation is valid. The rms value calculated for Cyg X-1 in the previous sections will 
then be correct within ^ 3% systematic uncertainties. 



